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Abstract 



c/3 i We present the complete formalism that describes scattering in graphene at low- 

■ ' energies. We begin by analyzing the real-space free Green's function matrix, and its 

analytical expansions at low-energy, carefully incorporating the discrete lattice structure, 
"j^ , and arbitrary forms of the atomic-orbital wave function. We then compute the real-space 

Green's function in the presence of an impurity. We express our results both in 2 x 2 and 

4x4 forms (for the two sublattices and the two inequivalent valleys of the first Brillouin 
^ . zone). We compare this with the 4x4 formalism proposed in Refs. [1, 2], and show 

O ! that the latter is incomplete. We describe how it can be adapted to accurately take into 

account the effects of inter- valley scattering on spatially- varying quantities such as the 

local density of states. 



oo : 1 Introduction 

Graphene's most interesting features are the presence of two sublattices and the existence of 
Q ■ linearly- dispersing quasiparticles close to the Dirac points, where the valence band and the 
00 , conduction band touch. The existence of two sublattices (or equivalently of two atoms per unit 
' cell) necessitates the introduction of an extra degree of freedom, often called pseudospin, which 
gives the Hamiltonian a 2 x 2 matrix structure. The Hamiltonian matrix can be diagonalized 
^ ■ to obtain the band structure, which is linear at low energies. 

j^ ■ This allows to extract analytically the low-energy physics by expanding the Hamiltonian 

close to the nodal points and focusing entirely on the regions in momentum space close to 
these points. However, restricting the analysis to low momenta disregards the real-space 
discrete structure of the lattice. For a monoatomic lattice, the discrete result can be recovered 
straightforwardly by overlapping the continuous result with the discrete lattice structure. 
However graphene's two sublattices make the overlap with the lattice less intuitive (a careful 
analysis of the real-space discrete structure of the Green's function is presently lacking). 

Here we derive the first complete real-space Green's function in monolayer graphene that 
takes into account the discrete structure of the lattice and the presence of the two sublattices. 
We find that different elements of the matrix Green's function are non-zero on different sublat- 
tice sites. We begin by analyzing the full real-space Green's function valid at arbitrary energy. 



This cannot be evaluated analytically, but can be obtained by performing a two-dimensional 
numerical integral. Our formalism allows us to describe both localized and extended atomic- 
orbital wave functions. Subsequently we derive analytically the discrete low-energy real-space 
Green's function, which we first write as a 2 x 2 matrix (using two sublattice indices). This 
is completely general and incorporates the existence of the two inequivalent valleys without 
requiring the use of 4 x 4 matrices (two sublattice and two valley indices). However, in order 
to make connection with the previously derived analytical 4x4 formalism [1, 2], we rewrite 
the Green's function in 4 x 4 language, and compare our result with the real-space 4x4 
Green's function derived in [2]. 

We then apply our formalism to describe impurity scattering. We use the Born approxima- 
tion to relate the free Green's function to the Green's function in the presence of an impurity, 
and to calculate the local density of states (LDOS). In the 2x2 formalism, the density of 
states can be obtained simply from the trace of the Green's function matrix. However, in the 
4x4 formalism the density of states is the sum of the traces of the four blocks that compose 
the 4x4 matrix Green's function. This modifies the 4x4 formalism of [1, 2], and allows it 
to describe correctly the effects of inter-nodal scattering. 

Indeed, for localized impurities (in agreement with [4, 5]), both the 2x2 and the corrected 
4x4 formalisms describe accurately the existence of the short- wavelength (-R3) oscillations 
[10, 11] generated by inter-nodal scattering. Moreover, these formalisms predict that these 
oscillations decay as 1/r [4], slower than the long- wavelength 1/r^ oscillations generated by 
intra-nodal scattering [1, 2], and in agreement with recent experimental observations [11]. On 
the other hand, the original form of the 4x4 formalism [1, 2], while capturing accurately 
the features of the intra-nodal Friedel oscillations (FO), cannot capture the short- wavelength 
inter-nodal FO for any type of impurity. 

In section 2 we compute the Green's function at arbitrary energy. In section 3 we analyti- 
cally evaluate it at low energy: in section 3.1 we approximate the form of the atomic orbitals 
by delta-functions and we compute the Green's function, writing it as a 2 x 2 matrix. In section 
3.2 we generalize this to finite-size orbitals. In section 3.3 we rewrite this Green's function 
in 4 X 4 language and compare it to the one proposed in [2]. In section 4 we use the Born 
approximation to compute the generalized Green's function in the presence of an impurity. 
In section 5 we calculate the dependence of the local density of states (LDOS) on position. 
We conclude in section 6. In the Appendix we present the derivation of the imaginary-time 
Green's function in momentum space. 

2 Green's function 

The tight-binding Hamiltonian for a graphene monolayer is given by: 

n = -tJ2ia]h + h.c.), (1) 

where t is the nearest-neighbor hopping amplitude, and (ij) denotes summing over the nearest 
neighbors. In momentum space this Hamiltonian can be written as: 

n= [ [aUk)Kik)fik)+h.c.] (2) 



where 

f(k) = -t{e-'^-^^ + e~'^-^^ + e-'^-^'') , (3) 

and the integral is performed over the first Brillouin zone (BZ). Also, 5i = av^x/2 + ay/2, 
82 = — av^x/2 + ay/2 and ^3 = —ay (as depicted in Fig. 1), where a is the distance between 
two nearest neighbors. As described in [3], there are multiple equivalent choices for f{k), 




^W^''^ 



a) 




Figure 1: Hexagonal honeycomb lattice of graphene (a), and its band structure (b). In b) the equal 
energy contours are drawn, and the Brillouin zone is indicated by dashed lines. The Dirac points K 
and K' are marked by arrows, and the reciprocal lattice vectors o^ 2 ^^^ ^'l^o drawn. 



depending on the tight-binding basis used, and on the definition of the Fourier transform 
(FT) of the real-space operators a^ and bi. The choice (3) (the // basis in Ref. [3]) makes the 
calculations more transparent and will be used throughout this paper. 

The real-space imaginary-time matrix Green's function can be defined as: 



Gap{f,iUn) 



^ {clj{iUn)Cpi{iuJn))5{Rf3j - R^i - f). 



(4) 



Ri,Rj 



where the sum is performed over all unit cell sites Rij, with Rj = nSi + ma2, and j = (n, m) 
specifying the position of one graphene unit cell. Here we take ai = av^x/2 + 3ay/2, and 

and 



-Ai 



lUJr. 



a] I lUJr. 



0,2 = — av^x/2 + 3ay/2. Also, a, P = A, B such that the operators 

c'^-{iujn) = bl{iuJn) denote creating an electron with energy iujn at sites of the sublattice A 
{Rai = Ri + Sa = Ri) and of the sublattice B {Rsi = Ri + Sb) respectively, where Sa = and 
Sb = S3 = -ya. 

The definition in Eq. (4) stems from assuming that the atomic orbitals are delta-function 
localized at the atomic lattice sites. If the atomic orbitals have a finite extent, Eq. (4) is no 
longer valid, and the delta-functions need to be replaced by continuous functions. This will 
be discussed in more detail in section 3.2. 

The corresponding momentum-space Green's function can be obtained by simply taking a 
FT of the real-space Green's function: 



Gap{k,iuJri 



(ffG{f,iuJn)e 



ik-r 



(5) 



The advantage of this generahzed Green's function is that it is defined for arbitrary values of 
k, and not only in the first BZ. 

We can easily show (see the Appendix for details) that this Green's function can be related 
to the two-point function of momentum-space operators: 

Gap{k,iuJn) = {cl{k,iun)ci3{k,iujn)). (6) 

While, by construction, the expectation value is defined only for k G BZ, as described in the 
Appendix, the above relation should be understood as evaluating the form of G{k) inside the 
BZ, and then expanding the validity of the functional form for arbitrary k. 

In order to obtain the retarded real-time unperturbed Green's function we calculate the 
expectation value with respect to the Hamiltonian (2) using an analytical continuation iun — *■ 
uj + i6: 

By construction, the real-space unperturbed Green's function can be related to the momentum- 
space G^{k,u) by: 

G°(f,c.) = |0G°(^,c.)e-*n (8) 

where we stress again that the integral is performed over the entire momentum space, and 
not over the first BZ. Dividing the momentum space in BZ's, we can rewrite the components 
of the Green's function matrix as: 



- „^ JkeBZ ^^ 

J2 I f^G%{k,uj)e-''^-'^e-'^^-^<'^+'-~'^^ (9) 



„, JkeBZ 

QfiueHL 



where al = x27r/v3a + y27r/3a, a2 = — x27r/v^a + y27r/3a, such that Q^^ = iia\ + va^ are 
the reciprocal lattice vectors, and /i, v are arbitrary integers. Here we used the fact that, as 
described in Ref. [3], G^^p{k + Q^y,uj) = G'^p{k,uj)e''^'"''^^"~^i^\ The sum over the reciprocal 
lattice vectors can be performed to yield 



d'k 



■^0 tu , ,\^—ik-r 



G^r, 00) = Y 6{f- R, + 6^- 6p) / -^G%ik, uj)e 

= J2 Sif-R,+l-6p)G%if,uj), (10) 

j={m,n) 

where G^{f, uj) is the reduced real-space Green's function obtained by integrating the momentum- 
space Green's function only over the first BZ. We note that the AA and BB components of the 
Green's function are non-zero only for the sites Rj of the A sublattice (we use the convention 
that an A atom sits at the origin of the coordinate system), while the AB and BA components 
are non-zero only for 5-type sites. 



Thus one can obtain the total Green's function by computing the reduced Green's function 
and multiplying it by the corresponding delta functions in the manner described above. This 
is quite general, and valid at arbitrary energy. However the integral (10) can be evaluated 
only numerically [4, 5, 6, 7, 8]. 

3 Low-energy analytical expressions 

3.1 The Green's function in 2 x 2 language 

At low energy, the Green's function can be obtained in a more elegant way by observing that 
only the regions close to the nodal points contribute to the FT integral: 



. 47r2 ^^' 



G\f,u;)= Yl I 1^2^10^ ^)e-'^'^''^^^-' (11) 

where [3]: 



7^ Ol On -,^ -^ 

K^V = ^ o ^ /^"l + '^^2- 

Here ^ = ± is the valley index (there are two such points for each unit cell of the reciprocal 
space), and r] = {fi, u), such that Q^^i, = fiSl + u(X^ G RL. For example, the most used valley 
pair containing two corners of the first BZ is described in this notation as -ft' = -ft'+(o,o) = 
(47r/3V3,0) and K' = i?-(o,o) = (-47r/3v^, 0). Also Gl^{q) = G%^{q + K^^), and g is a small 
deviation from a nodal point. Thus we obtain 

G\f,uj) = Y,G%{r,u:)e-^'^^^-\ (12) 

where 

Gl{r, c) = I 0G?,(g, c.)e-^^^n (13) 

The G^cAr^uj) functions can be calculated exactly by linearizing f{k) near the nodal points: 

where r = |r|, z = \uj\r/v^ v = 3ta/2, Hq {(z) are Hankel functions, and z^^^ = e"*^?'''^''')'*^'^-^-^^) _ 
g2j7r(^+i/)/3^ Noting that the sum over the reciprocal lattice vectors (/i, u) can be performed to 
obtain a sum of delta functions centered on the Rj sites, we obtain: 

G\f,u;) = YGlir,u;), (15) 

S=± 

with 



rO(^ , A ^ Y^ , ,^^ier-K ( ^{r- Rj)Hl Hz) i^r- Rj - '^s)(^a; - iy)H\ ^(z)/r 

(16) 



,=(„,„) V W- Rj + ^b){^x + iy)Hl'>iz)/r 5{r - R^)H'^'> {z) 



This can be evaluated to yield 

G»(r-„)c.2 V Js(r-Rj)H'^\.)cos(K.f, ^r- Rj - Ss)Mr)H[y \ 
J.M \'>(f-R,+SB)M'lnl^'(^) S{f-R,)co8(K-r)H'^\^) ) 

where 0i,2(O = [sin(_ft' ■ f)x ± cos{K ■ f)y]/r, K = -ftr+(o,o) = (^i ~ ^2)/3- The factors e*^'^ in 
the Green's function have also been discussed in Ref. [9]. 

3.2 Extended atomic orbitals 

The above derivation of the Green's function is valid if the atomic orbitals are fully localized. 
More generally, we can take into account the finite extension of the atomic orbitals by replacing 
the delta-functions in the definition (4) of the Green's function by a finite-extent function 0(r), 
centered at r = 0: 

Gap{r,iuJn) = ^ {clj{iuJn)cpi{iuJn))(l){R^j - Rai - r) ■ (18) 

Ri,Rj 

A possible choice for is a Gaussian (0(r) = e~^ /'«), where la is the "radius" of an atomic 
orbital. We can see from the derivation in Eq. (34) of the Appendix, that the FT of the 
real-space Green's function is no longer equal to a two-point function of momentum-space 
operators. For simple forms of such as the Gaussian, the integrals in Eq. (34) can still 
be performed analytically, but the two-point function in Eq. (36) should be multiplied by 
a momentum-dependent factor that decays at large momenta. In the low-energy limit, this 
translates into a decay of G°^(g, a;) with the magnitude of K^ri- Hence the nodal points do not 
contribute equally to the sum in Eq. (12), which cannot in general be evaluated analytically. 
Therefore the analytic form (16,17) of the Green's function can only be used when the extent 
of the atomic orbitals is negligible ^. 

If the extent of the atomic orbitals is finite, the Green's function can be obtained from 
Eq. (18). The two-point function of the Ca,i3 operators is simply the Fourier transform of the 
momentum-space two-point function (with the momentum integral performed over the first 
BZ), which needs to be overlapped with the atomic-orbital wave-functions using Eq. 18. 

However, this integral over the first BZ cannot be evaluated analytically. To obtain the 
low-energy analytical behavior of the Green's function one can start from Eq. (12) and use 
the fact that, for finite-size atomic orbitals, G°^(^, uj) decays with the magnitude of i?^^. This 
implies that the sum in Eq. (12) is dominated by the nodal points that are closest to the 
origin of the first BZ. Their exact number, and corresponding weights depend on the exact 
atomic-orbital function. The simplest approximation is to consider that only the nodal points 
at the corners of the first BZ contribute [4]. This yields: 

G'{f,co) = Y,Glir-co), (19) 



^If the orbitals are not fully localized, the Green's function in Eqs. (16,17) also breaks the hexagonal 
symmetry of the lattice. 



with 

''«<'^-"' ""1^ ,*,('') '^IK- + '!/)<«/'■ i^W'iX'W )• * * 

where #^(f) = (e-^^^-^e(o.o), e-^^^-^«(o.«), e-*^^-^«(-f.«)), 1= (1,1,1) and Z^ = (1, e^^'^^/^ e-^^^^^/^^^ 

3.3 The low-energy analytical form for the Green's function in the 
4x4 basis 

It is interesting to note that, instead of performing the sum over ^, one can also write the 
Green's function as a 4 x 4 matrix, with indices corresponding to the A_^, B^, A_, B_ compo- 
nents: 

G°^,(f,a;) = %G°(f,a;). (21) 

We can then compare the above Green's function (where G2(r*, cj) was derived in Eq. (16)), 
to the 4x4 Green's function proposed in Ref. [2]. The latter has a similar structure, with 
G^(r, w) being given by: 

G.I''--) «(iK, + ,:,)<(,)/, ffa)(,) j (22) 

The main differences between this Green's function and the Green's function we derived in 
(16) for localized orbitals are 1) the discrete structure, and 2) the presence of the oscillatory 
terms e^^'^ . These two factors are crucial to explain the short-scale fluctuations of the system. 

4 The generalized Green's function in the presence of 
an impurity 

To describe systems that are not translationally invariant (like in the presence of an impurity), 
one needs the "full" Green's function matrix, defined as 

Qaniri, r2, iUn) = ^ {clj{iUn)Ci3i{iUn))5{Rai - ri)5{Rpj - fs). (23) 

Ri,Rj 

For a translationally invariant system with one atom per unit cell, the relation between the 
usual Green's function and the "full" one is simply G^{f2—fi,iuJn) = ^°(^i, T2, iuJn)- Indeed, for 
a monoatomic system, the difference r*i — r*2 can be a Bravais lattice vector only if fi and r*2 are 
both Bravais lattice vectors. This yields^: '^- ■5{ri — r2 — Ri + Rj) = X]j j ^{^i~ Ri)^{'f^2 — Rj)- 
For graphene, things are however a bit more complicated. Thus, if 7^2 is on a site of type A, 
the components Q^^ are zero, while the components G^^ are not. In general the relation 
between Q^ and G" can be written as 



62^(n.r-i,i:u,„)= V G'ifAf2-n)i>{r,-B,„). (24) 



j=(m,n) 



^We use a normalization of the delta function in which S{r — Ri) = N ii r = Ri and otherwise, where N 
is the total number of unit cells in the crystal 



We would now like to use the above formalism to carefully calculate the corrections to the 
Green's functions in the presence of an impurity. In the Born approximation'^, the impurity 
corrections to the "full" Green's function are given by: 

Ag{fi,f2,iuJn) = / ^°(fi,f3,ic^„) ■^(fs) ■^°(f3,ri,ic^„), (25) 

•Jfi 

where ^(r*) is the scattering potential matrix, and the dot ('■') denotes matrix multiplication. 
At low energy, an analytical expression for ^° can be obtained directly from Eq. (17) and 
substituted in the above expression. However, it is sometimes more transparent to use the 
expansion (11) to write: 

Ag{n,f2, tUn) = Y. f ^".(^"l' ^"3' ^^n) ' ^(^"s) " Ql','{r3, ^2, ^O^Jc^^^^''-^^^-^^''''-^^^) (26) 

For r*! = r2 this equation allows us to see explicitly that each inter-nodal scattering process 
contributes at a wavevector corresponding to the difference between K^j and K^/j/. Also, we 
can easily modify Eq. (26) to approximate the spatially dependent matrix u{fs) by a matrix u 
which is independent of position, but that depends explicitly on C,,C',V and r]': ^(ra) -^ uf},^. 
In order to make connection with the 4x4 formalism described in [1, 2] we rewrite (26) as 

A^(fi, f2, iUn) = ^ A%/(fi, f2, iujn) (27) 

with 

A^55'(fi, fs, iUn) = / ^°(ri, fg, iUn) ■ Uin) ■ ^°,(f3, fs, tUn) (28) 

Jra 

where Q^ can be extracted from Eqs. (16,24). This corresponds to a 4 x 4 expression: 

A%'(fi, fa, iun) = jY ^°^j(fi, fg, iun) ■ %i52(r3) ■ ^^^/(fg, fs, iun), (29) 

where Gl^^ ^ 5^,^,01 (20). 

However, as can be seen from Eq. (27), the physical properties of graphene are not related 
directly to the 4x4 Green's function matrix, but to the sum of the four 2x2 blocks that 
compose it. Thus, the density of states is given by 

p{r, u) = } Im{Tr[^££'(r, r, iu;„ -^ u + iS)]}. (30) 

This differs from 



TX 



p{f, Lo) = V" Im{Tr[^^^(f, f, iu;„ ^ a; + i5)]}, 

TT -^ — ' 

which was proposed in [1, 2]. Together with the form of the full Green's function presented 
in Eq. (17), our prescription (30) is a crucial factor in describing correctly the effects of 
inter-nodal ^ -^ C scattering using a 4 x 4 formalism. 



■'Our analysis also applies to the T-matrix approximation if the impurity potential is localized. 

^For example, for a localized impurity, the scattering is independent of momentum, so u should be inde- 
pendent of CjC'i^ ^nd rj' . However, if the impurity is extended, the inter-nodal scattering is suppressed and 
■Uct, oc S^^'Sijrj' . 



5 The LDOS in the presence of a localized impurity 

We now focus on a simple limit when the impurity is localized on an A atom, (delta-function 
potential), for which we obtain: 

A^(fi, fa, iun) = ^°(fi, 0, iun) ■ u ■ ^°(0, r^, iun) (31) 

where « = {{!, 0}, {0, 0}}. For localized atomic orbitals this yields: 

p(f) = _iM[G%(-f, u;)GUr, uj) + G%{-f, uj)GUf, uj)] 

TV 

oc -i 5^ lm[6{f- RAj)Hl,'\zf cosiK ■ r)^ + 5if- RB,)M^^Hi'\zf] (32) 

j={m,n) 

Far away from the impurity (z ^ 1) the impurity gives rise to 1/r decaying oscillations on 
both sublattices. As mentioned in [4, 5, 12], the oscillations due to intra-nodal scattering have 
alternating signs on the A and the B sublattices. Thus the dominant 1/r behavior cancels if 
one coarse-grains the system so that the A and B atoms in a unit cell are measured together. 
This is what is expected in a real-space LDOS measurement in a real graphene sample, due 
to the fact that the A and the B atomic orbitals have a finite overlap. The A and the B 
oscillations are also mixed together when one measures the FT of the local density of states. 
Indeed, the intra-nodal oscillations give rise to a disc in the center of the BZ corresponding 
to 1/r^ decaying oscillations, and not to a ring-like singularity, as for 1/r oscillations [4]. 

For inter- nodal scattering, we can see that the 1/r decaying oscillations do not cancel 
between two neighboring A and B atoms [4, 5]. This gives rise to much more robust, short- 
wavelength oscillations, which can be observed for larger areas around the impurity. These 
translate into rings of high intensity in the FT of the LDOS close to the corners of the BZ. 
The difference between the intra-nodal and the inter-nodal FO was confirmed by a recent 
experiment [11], which observed no ring in monolayer graphene close to the center of the BZ, 
and detected rings close to the corners of the BZ. 

We should note that by using the formalism in Refs. [1, 2], the 1/r decaying short- 
wavelength oscillations are not accurately retrieved. For example, no 1/r decaying inter-nodal 
oscillations arise, even in the presence of a localized impurity. As we show here, this can be 
corrected by using the expression of the Green's function presented in Eq. (17), and the correct 
prescription (30) to extract the LDOS starting from the 4x4 "full" Green's function. 

6 Conclusions 

We calculated the complete low-energy free Green's function for monolayer graphene. We 
expanded carefully the Green's function at low energy, paying particular attention to incor- 
porate the effects of the discrete structure of the graphene bipartite lattice, and of the form of 
the atomic-orbital wave function. We have also calculated the "full" Green's function in the 
presence of impurity scattering. We have found that the 4x4 formalism proposed in [1, 2] 
to study the low-energy physics of graphene, while describing correctly intra-nodal scattering, 
does not accurately capture the effects of inter-nodal scattering on local quantities, such as 
the LDOS. We have shown that the missing pieces in the 4x4 formalism in [1, 2] are 1) 
the discrete structure 2) the presence of the oscillatory terms e"^^'^ in the Green's function 



and 3) the prescription to compute the physical properties not from the 4x4 matrix itself, 
but rather from the sum of the four 2x2 block that make up the 4x4 matrix. We believe 
it is important to check whether, besides the LDOS, there exist other spatially-dependent 
quantities that are inaccurately described by the incomplete 4x4 formalism. We have also 
presented a 2 x 2 formalism that is more general, and can be used straightforwardly for any 
type of disorder calculation, and for any type of atomic-orbital wave function. It would be 
interesting to use our formalism to calculate the real-space dependence of the LDOS in the 
presence of an extended impurity. 
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Here we show that the FT of the continuum real-space Green's function defined in Eq. (4) 
corresponds indeed to the two-point function of momentum space operators: 

Gaf3{r,iuJn) = "^ {clj{iuj.a)cpi{iuJn))5{Rpj - Rai - r) . (33) 

The FT of the real-space Green's function is 
Gai3{k,iuJn) = / d^fGai3{r,iuJn)e'^''' 

= d^r^ {clj{iuJn)ci3i{iuJn))S{Rpi - Raj - r)e'''"' 

Ri,Rj 

g g JkieBZ Jk2eBZ 

Noting that ^^^ e^^^"^')-^"' = Y.Q^,eRL^^'9{-^Q^lv■^a)^{k-k' + Q^,y), where Q^r, = fial + ua^ 
is any vector of the reciprocal lattice, 6a = 0, and 6b = (0, —ay), we have 

Gafs{k,iuJn) = {cl{k + Q^,y,iuJn)c^{k + Qf,u,i(^n))\k+Q^,eBZ ^ 

X exp[-iQf,i, ■ {6a - 6p)] . (35) 

We can see that this relation holds for arbitrary k, inside and outside the first BZ. Using the 
fact that cp{k + Q^u,i^n) = C/3(A;)e"*'5'"'''^'3, and cl^{k + Q^^jun) = Ca(A;)e*'3'"'''^° we obtain 

Gap{k,iujn) = {ci{k,iUn)cf3{k,iiJn)) ■ (36) 

We should note that if one uses the alternative (I) basis presented in Ref. [3], the different 
components of the matrix Green's function acquire different phases, consistent with the fact 
that the density of states written in this basis in momentum space is a superposition of 
two-point functions of A and B operators with a relative phase factor [3]. 

Another observation one should make is that if we repeat the above exercise for the "full" 
Green's function defined in Eq. (23), we find: 



= {ci{ki,iujn)cf3{k2,iujn)) ■ (37) 
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